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Sucima ry 

This final report covers work conducted by 3righaa Young 
University personnel in support of the Pioneer missions to 
Jupiter (10, 11), and Saturn (11) as well as the reduction, 
analysis and interpretation of magnetic field data obtained by 
the vector helium magnetometer (VHM) on the Pioneer 10 and 11 
spacecraft. Initially, our efforts concentrated primarily on the 
interplanetary data, and those aspects of the data of relevance 
to obtaining a better understanding of the interaction of the 
magnetized solar wind with the terrestrial magnetic field. 
However, after we participated in the encounters of Jupiter and 
Saturn, the emphasis of our research was directed primarily to an 
analysis of the planetary data. In particular, it soon became 
clear that there was a need for modelling of the various 
candidate magnetospheric currents suggested by the data. 

Over the period of this contract (7/1/73-12/31/81) we have 
supported the launch, cruise, and encouter phases of the Pioneer 
10 and 11 missions, published 16 papers, presented 15 papers, and 
had 19 abstracts of papers published. A listing of the papers 
and abstracts published is given in the "List of Publications." 

In the next sections we summarize the work conducted la various 
research areas conducted under this contract. Included are 
results not published as yet, but which we plan to publish. 

The August 1972 Solar Flare Event 

During August, 1972, considerable flaring activity occured 
on the sun, which in turrf produced significant modification of 
the solar wind medium as measured at a number of spacecraft. 
Including Pioneer 10. In addition, significant terrestrial 
effects were observed. We worked closely with E. J. Smith of JPL 
during this period and conducted some additional analysis of our 
own that were based upon the hypothesis that flare producing 
regions rather than specific flares are the causal source of 
Interplanetary and terrestrial events. Although we did not 
complete the analysis of this event in terns of the flare region 
hypothesis because of pressures resulting from the encounter of 
Jupiter, nevertheless there were a number of characteristics of 
the data (both field and energetic particle) that were consistent 
with this latter hypothesis. That is, effects were seen in space 
and at earth when the flare producing region (region on the 
surface of the sun having a history of flaring) was at the proper 
angular relationship with respect to the earth considering the 
finite travel time at nominal solar wind velocities, and it was 
not necessary to hypothesize that a flare erupted on the unseen 
portion of the sun. The results of a standard analysis and 
interpretation of these data has been published .although we are 
planning to develop and present our alternate hypothesis (L. 

Davis et al., 1 97 3 ; E. Smith et al . , 1 97 7 ). 



The MaasurenenC of Weak Magnetic Fields In Space 

During the period of this research we studied a proposed 
aethod to measure weak magnetic fields using two magnetometers in 
a coorrelative mode purported to reduce the resultant 
magnetometer sensor noise to an insignificant amount. We found 
that the basic premise of the paper was in serious error, and 
that in fact the method proposed would not work. The results of 
this brief study were published (Jones, et al., 1974). 

Interior Source Magnetic Field Modelling; Jupiter 

Our initial modelling efforts as applied to the Pioneer 10 
data consisted of the utilization of a function minimization 
approach to finding the field sources that best fit the data 
obtained within about 7 RJ (D. E. Jones and J. G. Melville, 1974; 
E. J. Smith, et al., 1974). When a single offset dipole was used 
as the interior source of the field, we found that as the radial 
range of the data used in the fitting was increased, equatorial 
projection of the end of the vector describing the offset of the 
equivalent dipole would tend to describe a circle. A more 
complicated field source configuration consisting of two dipoles 
plus a uniform field was then used, the latter being the dominant 
term of an external ring current. The use of two dipoles was 
based upon the results of ground based radio astronomy 
interferometer mapping and polarization measurements which 
suggested the existence of a magnetic anomally near the surface 
of the planet near a System III longitude of 200°. For 
simplicity, we chose to represent this anomaly in our studies as 
a magnetic dipole. The other dipole was assumed to be the 
primary source of the planet's magnetic field. 

A function minimization (FMFD) algorithm was developed which 
internally varied the location, orientation and strength of the 
two dipoles, as well as the parameters of the uniform field, 
until a best fit was obtained to the data in a least squares 
sense. Using the Pioneer 10 data, it was found that the 
resulting residuals were smaller than those obtained using the 
conventional linear, or spherical harmonic, approach. We found 
that a plot of the residuals versus radial distance exhibited 
high values at about the orbital distance of lo . The Voyager 
later confirmed that the flux tube of lo contained currents, and 
suggested the presence of a plasma torus at this orbital 
distance. The second dipole was found Co be located close to the 
equator in Che hemisphere, close Co the longitude of the radio 
astronomy polarization anomaly, and reasonably close to the 
surface. 

A similar type of analysis was performed using the Pioneer 
11 data with results Chat were consistent with that obtained from 
Pioneer 10. The second dipole was located in the same 
longitudinal quadrant as Chat of Pioneer 10 and reasonably close 
to the surface, but in Che southern hemisphere. We are not sure 
what effect the existence of tail and frontside cross sheet 
current fields will have on the location of Che second dipole as 
determined from the Pioneer 10 and 11 data, but it is clear chat 
Che use of such an analysis approach has merit, particularly if 
proper allowance for the fields from these other currents results 
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in no re consistent characteristics for the two interior dipole 
sources as derived froti the two spacecraft data sets. An 
extension of this work, which included the Pioneer 11 results as 
well as a SHA analysis of data from both spacecraft, used a 
different weighting scheme, and derived more information 
concerning the ring current, was written up as an undergraduate 
thesis for the honors program by R. Steven Turley and has been 
included as Appendix A of this report. 

Modelling the Jovian Magnetosphere 

During the period of this contract we have studied the 
Jovian magnetosphere using two basic methods. The first involved 
the use of Euler potentials, a method which results in 
mathematical expressions that permit easy tracing of field 
lines. Some degree of success was achieved in obtaining 
functions describing the field which provided reasonable fits to 
the data, and which were useful when extrapolated a small 
distance beyond the region of fitting. As the result of a 
cooperative program with the energetic particle experiment team 
at the University of Iowa, it was found that the outbound Pioneer 
10 was periodically located on open field lines, this occurring 
while at fairly low magnetic latitudes (Goertz, et al., 1976). 
This result was obtained by comparing the spacecraft location in 
magnetic coordinates with the energetic particle measurements. 

The existence of an open/closed field line demarcation at such a 
low magnetic latitude (approximately 20°) was rather startling at 
the time, but this may be consistent with the general topology of 
the magnetic field in the sunward magnetosphere as inferred from 
the outbound Pioneer 11 magnetic field measurements, which 
suggests that it is quite different from that of earth. 

We attempted to use an Euler potential approach in the study 
of the magnetosphere beyond ± 20° ma gne tic latitude, with little 
success. We found that the unphysical magnetic fields and 
currents predicted at high latitudes by the functions that were 
tried far outweighed the advantages of these functions (Jones and 
Melville, 1975). In addition, the deformations of the current 
disc evident in the data could not be accomodated in a tractable 
manner with such a function. We tried several coordinate 
systems, with little success. 

In order to facilitate a current disc clearly displaying 
twisting and bending, we subsequently developed an algorithm for 
the magnetic field of a double layer of circular rings of current 
which allowed the tilt and longitude of the axis of each ring 
making up the disc to be functions of the ring radius. In 
addition, we desired a model for the currents which could be used 
to extrapolate beyond the region of fitting with more reliability 
than the Euler function method. Excellent agreement was obtained 
between model predictions and the data, with parametric 
expressions being developed for the manner in which the current 
disc was deformed, particularly using Pioneer 10 outbound data. 

It was found that the disc needed to be twisted about the spin 
axis, and bent such as to approach parallelism with the 
rotational equator. An additional deformation of the current 
disc in the form of a small ridge or hump was required in order 



for the spacecraft to periodically penetrate the bent disc. 
Figures 1-6 compare the current discs required to fit the Pioneer 
10 outbound data (Figs. 1-3; twisted and bent, with spiralling 
ridge) and the Voyager data (Figs. 4-6; twisted and bent only). 
The Pioneer 11 inbound data also suggested that some degree of 
bending improved the fits. The other data sets clearly required 
additional tail field and required discs displaying markedly 
different characteristics as to the amount and distribution of 
current, as well as the amount of tilting required (Melville et 
al., 1975; Jones et al., 1975; Jones and Melville, 1975; Jones et 
al., 1976a, b; 1980; 1981; Jones and Thomas, 1931). An extreme 
case in point is the fact that in order to fit the strong current 
dip observed outbound at about 8 Rj by Pioneer 11, a tilt of over 
35*^ is required, whereas the inbound Pioneer 11 data require a 
thinner disc having a tilt of 11.2°. A study was also conducted 
to determine the radial current distribution required to produce 
a fit to the azimuthal field from both the outbound Pioneer 10 
and 11 data. 

When the characteristics of the ring currents for the 
various data segments we re compared, both the rho and phi field 
components displayed a local time dependence. As an alternate to 
local time dependent disc and radial currents, it was suggested 
the possibiliy of an azimuthally symmetric disc current plus a 
sheet of dusk to dawn current, the latter being the equatorial 
portion of a tail-like current configuration extending into the 
frontside magnetosphere (Jones, et al., 1981; Jones and Thomas, 
1981). Such a configuration produces the observed maximum in the 
rho and phi components near the dawn line, and the decrease in 
these components near the noon meridian. A preliminary study of 
the fields of such a current configuration has resulted in fits 
to the data which are better than those obtained in terms of a 
local time dependent disc current and t ai 1 /mag ne to pause currents 
(Jones and Thomas, 1981). 

Magnetic Field and Energetic Particle Flux Studies 

Several analysis efforts were conducted in an attempt to 
better understand the relationship between the energetic particle 
flux and the magnetic field. Using the model disc current which 
best fit the Pioneer 10 outbound data, we traced out field lines 
to determine the L values corresponding to important times 
related to observed characteristics of the energetic particles. 

As a result, a joint paper was presented outlining the L 
dependence of the energetic particles in Jupiter's magnetosphere 
(Jones and Mihalov, 1978). This study was based upon a model 
disc current system fitting the outbound Pioneer 10 data, where 
good fits were obtained without tall and magnetopause currents 
(these currents produce primarily positive and negative z fields, 
respectively, and they are apparently of the sane magnitude near 
the dawn meridian). Hence, the magnetosphere was not terminated, 
and a reanalysis of tiie particle data using L values derived rrom 
a terminated magnetosphere, or using model paramters derived rrom 
fitting other data segments should modify the results reported 
previously. 

We have also explored the locations of minimum field reti 
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that result when a magnetosphere is terminated by the solar 
wind. A model fitting the Pioneer 10 data has been combined with 
a system of currents on the surface of a sphere to derive a 
terminated magnetosphere. The resulting magnetopause was blunt, 
and consistent with the derived shape of the magnetopause using 
minimum variance techniques. The field minima produced by the 
termination were shifted away from the magnetic equator, 
resulting in two surfaces oriented at a relatively large angle 
relative to the magnetic equator. When the Pioneer 11 outbound 
trajectory was superimposed on the field line disgram derived 
from the model, it was found that the occurances of the field 
minima were reasonably consistent with the observed occurances of 
the energetic particle maxima (Jones, 1979). Other analyses of 
the Pioneer 11 outbound field data suggest several other possible 
explanations for the anomalously high flux counts at such a high 
latitude. One considers the possibility that the localized tilt 
of the current disc near the noon meridian was much greater (> 
30*^) (Jones et ai., 1975; Jones et al., 1976a, b) while the other 
considers the possibility that the magnetospheric cleft is at a 
much lower latitude, the latter providing a region for the 
energetic particles to escape from the planet (Jones and Thomas, 
1931) . 

Satellite-Magnetosphere Interaction Studies at Saturn. 

When the 24 hour interval of data spanning the closest 
approach of Titan's orbit by Pioneer 11 was studied, it was found 
that there were three periods in which the level of magnetic 
variability was enhanced^ One of these occured spanning several 
hours before to one or two hours after closest approach to Titan, 
when the spacecraft was about 145 (Titan radii) ahead of the 
satellite in the direction of orbital motion. The 
characteristics of these turbulent regions were studied and it 
was concluded that the Titan interval displayed a number of 
characteristics that are consistent with what would be expected 
should the spacecraft penetrate the satellites magnetic wake. 
Evidence for penetration of a shock so far "downstream" was weak, 
and marginal at best, but there was a field minimum almost 
precisely at closest approach to the extended tall axis, and the 
characteristics of the magnetic turbulence of tne Titan Interval 
appeared to differ from those of the other two. As a result of 
this analysis, we suggested Chat the magnetic wake of Titan may 
have been detected (Jones et al., 1979; Jones et al., 1980). 

Preliminary studies of Che magnetic data near the time of 
Che energetic particle decreases attributed to the new satellite 
1979 3-2 suggested the possible existence of low Alfven mach 
number fans associated with the interaction of a satellite with 
corotating plasma. However, a closer inspection of the manner in 
which Che various field components varied during this interval 
suggested this not to be the case. It was later determined that 
this was a data anomaly chat resulted from the use of a despin 
function that did not properly take into account acceleration of 
the spacecraft at Saturn and the nonlinear characteristics of an 
improper despin function timing anomaly in Che data produced the 
false result (Jones et al., 1979a, b). 
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Figure Captions 


Figures 1-3. Three views of Che currenc disc required to fit 
the Pioneer 10 outbound data. Twisting, bending (Cowards 
parallelism with the jovigraphic equator) and a deformation 
in Che form of a ridge having a spiralling symmetry line are 
required. The three views represented in the figures differ 
by 45° in the direction of rotation. 

Figures 4-6. Three views 45° apart in the direction of 
rotation of the currenc disc required to fit the Voyager 
magnetic field data based upon fitting of the outbound 
Pioneer 10 data. An additional deformation in the form of a 
sprialling ridge was not needed for the Voyager data. 
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Two models have been used to fit Pioneer 10 and 11 magnetic field data 
from Jupiter. The first model consists of two offset dipoles with an exterior 
current ring or disc. The second model is based on a truncated spherical 
harmonic expansion of the magnetic scalar potential v,'ithin the region of in- 


terest. Data used in the fits was obtained between l.£ and 7.0 Rj, 32° and 
-50° latitude, and 0° and 360° longitude. Data for Pioneer 10 and 11 were 


run separately and together. In cerms of an offset dipole and ring, the field 
is best described as a dipole with moment 4,1 Gauss-R,^ tilted 9.5° with respect 

U 

to the spin axis and offset about .14 Rj from the planet's center. The exterior 
current disc has an axis roughly coinciding with that of tne main dipole and 
an inner radius of less than 9 R,. 
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i^uRODucTio:; 

Ws nave fitted two models to Jupiter's magnetic field based on Pioneer 
10 and 11 data between 1.8 and 7.0 R,. The first model is that of two offset, 
tilted dipoles, with an exterior current ring. The second model is a truncated 
spherical harmonic expansion of the magnetic scalar potential. Data for Pioneer 
10 and 11 have been fit both separately and together in the cases of both 
models. 

The first model -was chosen because of how easy it is to visualize physically. 
Owens Valley Radio Astroncm.y data (Eerge and Oulkis, 1S76) and measurements by 
others including Branson (1963), Warwick (1964) and Conway and Stannard (1972, 
1976) suggest that the interior field sources are a main dipole, tilted about 
10° with respect to the spin axis, and a possible anamoly at longitude 200°. 

The ring current is ah attempt to model the current disc which extends 
from about 10 to 80 Rj and its contribution to the magnetic field near the 
planet. It is acmitedly not as good a model as, say, a current disc would be, 
but the contribution to the field here is so small that the ring model does a 
sufficiently good job. 

There is justification for including all three ccmpcnents in the physical 
model. The first justification is the evidence from other sources mat tnere 
are at least three contributions to the field in this region. The second 
justification is that there are obvious anamolies in the fit if any of these 
contributions are left out. Thirdly, addition of any of these three contri- 
butions to a model consisting of just the other two s igni fi cantly improves 
the goodness of tne fit. 


1 he 

second model 

fit 

ted to the field data 

is a truncated sphe>‘'cal 

harmoni c 

expansion of 

the 

scalar potential fro 

m whicn tne magnetic field 's 

ODta '! nec. 

This model 

has 

the disadvantage of 

being a little harder to 
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visualize p.nysically. las advantages are tnat it is a more general aDoroacn, 
and that it is a linear formula, vvfiicn makes finding a least souares fit con- 
siderably easier. The details of this approach are explained in more detail 
in a later section. 

In the cases of both models, we were able to obtain fits to the data 
with from .3% to 1% deviation from the experimental fields. We found the 
fitted values to be quite sensitive to the segment of the trajectories used 
in the fit end to the weighting factors which we used. We believe the 
particular weighting function used for data reported in this paper to be tne 
most reliable one. Interpretation of the differences between the various fits 
will be Given in a later section. 
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The experimental magnetic field data was gathered by a vector heli;.- 
magnetometer aboard the Pioneer 10 and 11 spacecraft. A detailed discussion 
of the instrument and its associated errors can be found in Smith (1975). In 
this section we briefly describe the experimental uncertainty in the measure- 
ments obtained and its relationship to the value of tnat is to be expected 
for eacn measurement range. 

Tne magnetometer has an absolute uncertainty of less than 5;i based on 
inflight calibration data. For the large fields measured in the experiment 
however, resolutions of at least ±lv were possible (.Ols: uncertainty). In 
this case the primary source of error is due to digitization. The following 
table was taken from information supplied by Smith (1974) and gives an indica- 
tion of the variance expected in each of the magnetometer's field ranges due 
to quantization of the data. 


Maximum Field 

Digital Step Size 

y/bit 

C" (■■■} 

4 V 

0.015 

1.88x10"- 

13 V 

0.052 

2.25x10"" 

43 V 

0.167 

2.32x10"- 

146 V 

0.569 

.027 

632 V 

2.46 

.504 

3380 V 

15.1 

19.0 

0.227 G 

83. 2 

648 

1 .37 G 

531 .0 

^ 3b 000 


is found as follows. Let e be the diaital step size and B be the 

0 

digital value of the field, B. If B is between - c/2 and -r c/2, the 
reported value will be 3. Assume equal probability of B being anv^where in tne 
interval. Let x be the difference between 5 and B,. 


— "1 ' - 

- ; / 1 . 






Note that on the average, a is of the crcer of .2:: of tr.e "ici-coint 
value for each range. This should represent the lowest possible value o 
c obtainable in any model fitting to the data. 



This section describes the units used in the calculations and results. 

First the position coordinates then the magnetic field variables will be dis- 
cussed. 

Throughout this oaper all position variables will be expressed in itodified 
System III (epoch 1957.0) coordinates, unless noted otnerv.’ise. Soherical 
coordinates are expressed in terms of radius, latituae, and longitude. The 
radius is in meters or in R, (Jovian Radii). Throughout this oaoer the value 
of 1 Rj = 7.08x10^ meters has been assumed. The latitude measurements are 
identical to System III latitudes. The longitudes are measured from the same 
axis as in System III, but in a counter-clockwise direction as viewed from the 
North Pole rather than in a clockwise direction. The modified System III 
longitudes (e') used here are related to System III longitudes (-jjj) by the 
relation 

e- = 360° -6jp. 

Dipole offsets are expressed in cartesian coordinates, using the spin axis as 
tne 2 axis and the line of 0° latitude and longitude as the x axis. 

The magnetic field variables are needed to express the dipole strength, the 
field strength, and the current in the ring of Xodel III. The dipole strength 
is exoressed in Gauss -R,’, with R, defined as before 

u 0 

1 Gauss-Rj- = 3.55 x 10^° Gauss-m’ 

= 3.55 X 10*° Weber-m 

Magnetic field data is in y, v/here ly = 10”- Gauss. Progrem calculations and 
outputs are in Gauss (1 Gauss = 10~~ Weter/m-). The currents are measured in emps. 

It will be noted tnat our choice of a fixed olanetary coordinate syste.m 
differs from that of Smith (1974), in which a coordinate syst=~ fixec to tne 
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Jupirer-Sun line for the external sources was chosen. The outer eoges of rne 
external sources is probably very sirongly influenced by tne solar wind, 
however, the dominant force shaping the near exterior field sources (whicn v«re 
of primary interest in this study) should be the strong planetary field. A 
partial justification of this assumotion is in the improved fits ve viere able 
to obtain compared to those reported by Smith et al , (1976). 


/ 



Nu';-Lir;£AR yoPEL c;lcI)LATIc:;$ 

In this section the formulas used to calculate the fields due to an 
offset dipole, a current ring, and a current disc vere explained. The 
formula for an offset dipole is vvell known and hence not derived. The 
formula for a current loop and a current disc are derived in detail. 

Dipole Field 

The equation for the magnetic field due *o a dipole centered at the 
origin is 

t = u^/4T:[M/i7l3 - ({M.?)?)/!r|S] 

where 

= 4ttx 10'^ (Gaussian units) 

M is in Gauss-Rj^ 
r is in Rj 

letting C be the offset of the dipole and = ^ - t, the above equation can 
be expressed for an offset dipole as 

B = u^/4-[R/';^;2 _ 5] 

E,^= y^/4-[M^/l^i3 - ((il.^)RCj/;^i2] 

B = y /4 -tLM /:^!- - ((M-^)RC ,)/;^! *J 
Z 0 z ■ Z ’ 




It is assLi~ed tnct r,?.s do J coTOonen", sines s^.'nrstry v;culd sesr^ to 
reauirs that 6 not depend on tne : coordinate of x. 

The 3 iot-Scvart Law cives o(?') = 7^ I C ^ 

" ' ht: 1 I -p- . 


I I - 

! To - i - 


i"2 = ri + 2k fi = a cos ©i + a sin ©j 

fs “ rj = (r - a cos c)i - a sin “ij + zk 

= {r-- 2 ar cos : -i- a^ - 
= (a^ + -r - 2 ar cos - 

dz = -a sin c d:i + a cos - dej 

dX X (r2-f.J = d? {a sin^ ik + z sin ij + (a cos- : - r cos i) k + z cos 

= ad-{z cos ?i + zsinij-^(a-r cos ;■) k } 


b = 


t la ^2r 
0 


I z cos C-; -i- z sin rj T fa-r cos :)k 

4- 1 (at + 4. 2^ . 2ar cos i)-/^ 


Since the denominator is even and the d 3 is odd, but d 3 and dE are 

y X z 

even, this can be written as 


u.ia r- 


B = 




f" z cos bi -i- (a - r cos k 

(at + + z- - 2ar cosn:t/t 


d: 


Lets = -JZ . 


la rT./Z 


B = 


z cos 2 =i + ( a-r cos 2 r' k 


df 


(at + rt + zt - 2ar cos lb) 




B = 


cos 26 = T - 2 sint e 


“0*^ (z-2z Sint x (g.p X 2r sint 

' •' (at rt + zt - 2ar - 4ar sint e)t/t 


de 


[">' ^ (z - 2z sin- 6)i - (a - r 
I ( 1 ^ 1 - sipt e)t/i 


zr sin- 6 ,K 


C6 


wnere 


(a- 


~ar 
-'i z 4 . 


, 3-r ' - 



rhi£ involves integrals of r.vc forr'.s 


A{£) = 


. r 


’ll 


d5 


J ( l+i-^sin^e) 
0 


and B(£) = 


:t/2 


sin'^ ed: 


J (1 + sin^9)“/^ 

0 


B = zG fA(£) - 2B(£)j i + [G(a-r)A(£) + 2rGB(£)] k 
To evaluate A(£) and B(t) let t=sin £ cos 5 = /I -t^ 


.1 


A( 0 ~ 


at = cos sae 
dt 

' , / r 


G£ = 


dt 


I (l-t^)-/^(l+£‘^t^)“/^ 


B{£) = 


V' i - 1 ^ 
j 


4 . 

U G « 


I (1-t^) V^(l+£^t^) V: 


1 

/•* 


dt 


A(0 = I (tl+i-2)[(t2+i-2)(l-t2)j- 




if a = Y b = 1 X = 1 This corresponds to 17.4.51 in Abranowitz 
and Stegen (p. 596) . 

1 


sin- = 


£^ 

1(t2+1) 


jr_ 

2 


k-^ = 




£- + 1 


_ 1 r- f r ! I 1 w 1 . _ [Ef £“/£-+! )j 

A(0 - T3 E i Ti 7T^:L( 72)( Y2 - ‘ .. ^ 

Z , A, Tl^, X, A, V-£-+l 


3(t) = r .3 


t^dt 


= 1. r dt 

I (t2+£‘2)[(t2+£’2)(]-t2)]l/2 j [ ( 1 2+^- 2 ) ( 1 - t2 ) J 1 / 2 


/ A ' 

i-A -^1 


With a, b, and x as before, this also corresponds to 17.4.51 in Abramov/isz 
and Stegen (p. 596) . 

2 1 


B(i) = T3 f 


2 |[ T2 + ^ 


I 


A "V A ■ 1 


A ■ 1 A- ■ I 
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Virere E(i) 


(V 


• dt 


is the complete elliptic integral of the second kind and 




K(0 = C(l-t2)(l-U2)] 


. 1/2 


dt 


is the comolete elliotic inteoral of the first kind. 


The Current Disc Field 


The equation for the field due to a disc are found from those of the ring 
in a simple manner. The equation of the ring current can be written as 

B = IA(r,5,i)) 

, , .T J dr 
let al = _o 

ri . 3 



where r is the radius of the ring 

**0 is the inner radius of the disc 
r^ is the outer radius of the disc 

The above integration was performed numerically, using the form.ula for 
J. derived from the formula for B derived in the previous section. 
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SPHERICAL HARMONICS CALCOLATIGNS 


If there are negligible currents in a region of interest, and if the 
electric field is slowly varying, v x B '\- 0, and B can be written as the gradien 
of a scalar. Let U represent the magnetic scalar potential defined by B = -vU. 
Since v • B = 0 requires tnat = 0, U must therefore satisfy Laplace's 
equati on. 

A solution to Laplace's equation in spherical polar coordinates is the 


familiar spherical harmonic expansion 
" A 
; \ } 

t=lm=0 


^ I Z Lr"''"'*(g'^ cos m i + hZ sin m cos m c + hZ sin m ;)] 


(cos 5), 


( 1 ) 


Here P'l (cos 9) are the Schmidt normalized Associated Legendre Polynomials 

I, 

defined bv 


P^(x) 


2(t-m)l 


( £+m ) ! 


V2 ^m 


Consider a region bounded by two concentric spheres. The contribution to 
the potential inside this region can be due to sources interior to the smaller 
sphere (interior sources), or exterior to the larger sohere (exterior sources). 
The coefficients gZ and R'Z in the expansion of the ootential due to interior 
sources must be zero if the potential is to be finite as r aporoaches infinity. 

m rn 

Likewise, the coefficients g^' and h^ must be zero if the potential due to 
exterior sources is to remain finite at r=0. For this reason, c'" and hZ' will 
be refered to as interior coefficients and cZ! and h'!' will be refered to as 

_e e 

exterior coefficients. The terms with t=l are comnonly called dipole terms, 
those with i=2 are called quadrupole terms, and those wish •'.=3 are refered to 
as Gctupole serms. 
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SDnerical Harrionics 


Many of the coefficients in spherical hanr.onic expansion can be related 
to parameters of our other model. In the following sections we will show how 
the coefficients are related to the parameters of a single offset dipole and of 
a current disc. 


Offset Dioole 


The potential due to an offset dipole of moment m, at a di solacemenc r 
from the dipole is: U = • "*st M will be in Gauss-Rj-). 

M • j" 

U = ^ — . Consider a coordinate system with an origin 0 which is different 
r^ ^ 

from tne center of the dipole, C'. Let r be the vector from 0" to the point 
P where the potential is to be computed. 



r = r - r' 
0 


Le* 




f = r sin 8 cos :i + r sin 8 sin oj r cos 8k 


Therefore , 

M (rsine cos:-r,) - M irsir.8 sin8-;' + *i_(rcos8-8) 

u = ^ ^ TT- 

[(rsine cosr-',)- (rsine sin:-:)^ (rccs8-8 
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The linear independence of the tems in the sonerical nantonic exoansicn 

require that coefficients of equal powers of r in both expressions of U 

(Equations 1 and 2) be equal. To enable this comparison, expand the denominator 

of equation 2 by the binomial theorem. 

After a little reduction, r can be written as 

c 

r- - 2r(^cose + nsine ccsc * c s inesin<j )]"-/" 

letting b = r,^ + q- - 

a = -2('COse + n sinecos; + ; sinasinr) 

+ ar + b] = [r^ + cj where c = ar+b 

^5/8c^ . 

r" • • • 

15/g,^ . 3/2b ^ 

• • • 

r5 

Since we are only interested in dipole and quadrupole terms in the 

expansion, only terms up to order Vr‘* will be retained. Hence, 

'13 - 

U = ^ • — :r(?cose + r,sin6cosc + ^sinesin;-)} ;M (rsinecosc-ri) 

r"" r ' X 

My(rsinesin(i-q) + rcos5-t ) } 

Ecuatinc -A— terms cives 
r^ 

M sinecos:: + M sinesini + M cose = a'xose+G‘ccstsine-^h^sin:sine 
X y '1 ''1 

Therefore , 

give 

r-M C* 3 (cCOSr-^nsi necos;--Si nesinc ) (>', sineccs:-’: sinesin 

X V z ' ' ' ■ ' X y 



M =0' 

X ■'i 


= h^ 


H 

V 1 

A 

^ = 5? 

- I 

The terms in 


[r- + c1‘2/2 = -3 - 

1- - y.0 


3/2c 


3 /o , 
/ ea 


U 



g- (3 cos^e-'i); - /T q-sinecosrCGS:- v''5 

./*T '’T 

' ^ q^siri'scos 2c ^ h-sin^rsin 

This equality must hold for all e and c; therefore it must be valid when 6=0. 

In this case 

q= = -r,H -;M +2cM, 

~2 X y ^ 2 

or, in terms of g^, g- , and h^ 

q-' = -rg-^-^h'mZ^q “ . 


n*STncST necosr 
2 c 


The coefficients of sinecosecosc and of sinaccsasine must also be equal, 
if the above exoression is to be valid for all 8 and This yields, 

gl = /J rM + /3n?^, 

2X2 

or 

g’ = /T(cg^ + ng°) 

^ 11 

and 


hi = /3rM^ + 

= /3(chi + cg^) 

1 1 

Substituting these expressions into equation (3) yields 

' q^sin^3cos2'b+’^ ^ h2sin^3sin2c = !' sin-3- \ ;M sin 

2 *2 2 2 2r, X - j 


3H-3y sin^scos-c 

XV 


3M tsin^ssin 

y 

Solvino for and h^ 
"2 2 


+ 3M nsin 

V 

w 

yields 


cosisinc + 3My ;sin-3coscsinc 


g- = *''3(n>! -?M^) 

1 ^ y 

h- = /3(r,iV| +cM ) 

2 y X 

or g- = /3(ng^-th* ) 


'3(r,h-+2g^ ) 


1 r 


CO 



Following the derivation by Schmidt (1934) v;e will find tne values of 

M and which would best represent the spherical harmonics expansion of a 

dipole field. Note that if the dipole were in a coordinate system with an 

origin at O', 5=n=4=0» and all and h”)) with t>l would be zero. To best 

approximate an offset dipole field, we will require that in the coordinate 

system centered at 0' [(Uo)^! be a minimum, where 
2 

U 2 = cos mi + h*^ sin m 0 ) f”)' (cos 5) 

[iUzPi = T I (s”) + (h") 

^ m=0 - 2 


minimizing |(U 2 )^| requires that vj(U 2 )^| = 0. 
In tne new coordinate system 


g^' = 0° 

1 “1 


g^' = g^ 


hi' = h^ 
1 1 


gO- = gi - (2cg0 - ng^ - thi ) 
22 111 

gl' = gl - /J(^gl + ng°) 

2 2 11 


hi' = hi - /3’('hi + cg°) 
22 11 


g2' = g2--/3(r,gi - thi) 

22 11 

h^' = h^ - v'?(rih‘‘ + tgi ) 


5i(U2)‘i = 5j(L)2)^' - 2'(2g2g0 + /Jqigi4>/3'hihi ) 

12 "12 12 

-2n(-gig^ + /Sg^gi + /3gig2 + /Thih^) - 2t(-hiqi -r/3n'hi 
12 12 12. 12 12 12 


- /Shiq^ + / 3ai h 
1^2 " 1 


+ (2^g^ - ngi - thi + 3[(?gi + ngi)' 
111 1 : 


+ (nhi + cgi)"J. 


(;hi + -i- (ngi - thi)i 

11 11 


I et 


A = 2a “O'-' -r v'3(g'‘g- h*h",; 

S = -G^a" + v-'3(0'io‘ -f- nlQt -i- ^lf-|2 

■'rt "r: ■'I't 1 2 

C = -hia" T/3(cth^ - 0 ^ 0 ^ + aih^) 
1'2 "1 2 i'2 '1 2 


} 



SinDlifying the aoove expression 'leads to 

5i(U2)-| = 5i(U2)^i - 25A - 2nB - 2'C - ^ r 

Setting the gradient of the above expression to zero requires that 

A = + g° (?a° + nd^ + th^ ) 

1 ■ 1 1 1 

E = 3M“ri + g* + r|G^ + th^) 

1 '1 "l 1 

C = T h^ (?9” ■*■ nO"^ + th"). 


(4) 

(5) 

( 6 ) 


Si nc£ 


g=A+g-B^h-C 
1 1 1 


3M2(5gC+r,g2+-h^ ) 


(c^'2+0^^+h^2) (rg^-rro^-rrh^ ) = 4M^ ( - 0 ^+r|C - + S h ‘ } 
"1 : i *1 ' 1 '1 "i 


then 


+ ”0^ + ?h^ = 

"1 ’1 1 


g?A + gjB -i- h’jC 


= D. 


With these definitions, C,n, and t can be solved from equations (4), (5), and (6), 


? = 


A - g?D 


3iw 


B - g^D 


; = 


C - hiD 

~~W 


Disc 

The spnerical harmonic expansion of a disc can be computed from the 
expression for the potential on the axis of a current ring. Let z be the 
distance along the ring axis, and a be the radius cf the ring. It is easy 
to show that 


U = 1/2 ll - 


a^+z‘ 


Let I have a radial dependence, dl = 


a- da 
c 0 


dU = 


J a“da 
0 0 

23^^ 




The bracketed terms can be expanded to give 
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~ J 3 
0 C 


_ f. 


-J 


Z J a ■ 

C 0 


If the axis of the disc is parallel to a spherical-polar axis, the 
spherical harmonic expansion becomes simply 


U = y r^A P, (cos e ) 
tSl ^ ^ 


at 6 = 0 


cc 

U = T 
t=l 


The coefficients of Z must be equal, therefore 


Ao = 


a 

0 0 


A, = -(^)J 


0 


Ao = 0 


A3 = (^) 


To transform the coordinates into a system not with the ring, the Addition 
Theorem is employed. Using Schmi dt-normal ized Legendre Polynomials gives the 
simple fort 

P (cos e") = I pT(cos £) p'['cos(ti) cos m (£-5) 

^ m=0 

where 5 ' is the angle made with the ring axis, and n and ; are the angular 
coordinates of the ring axis in the new coordinate system. Thus 

[cos cos It; - sin mt sin m;j 


y r"'A.F‘r (cos 6) F"' (cos n) 
£=1 m=0 " ^ 

r 

1 COS m.- ^ h.' 

m=0 
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iri 


ar,d A,P 7 (cos r,) cos mi = c. 
A,P''‘(cos n) sin mi = h 


X 


Comparing coefficients we obtain 


gy 


9 = 


Gt = 


/1+Olv « 

(— ) cos n 

V— 5— ) Sinn cos.; 


^5 = 

93 ° = 


- (^) Jq Sinn sini 
“^0 " 0 " 


3 cos n)) 


Solving for J , a , n, and 5 yields 
o c 


-1 

i = tan ^ ■=■,- 

0 T 


/hJ^+gp 

n = tan‘i [ i__ ] 


?o 

v» ^ 


0 - = 


0 1+ct /gp+gp+np 

_ ^ . (3+a)Jo(5cosp -Scosn) : -j ^ ^ ‘ 
8g? 


a^ will be in R, 

0 J 

J will be in Gauss 
0 


Conclusion 


In conclusion, in a current free region with stationary fields, the 
Doteniial can be expressed as a linear ccnbinasicn of soherical harmonic 
functions. Using the above formulas the coe'fi cienos of onese functions can 
be related to the model oarameters of a single c-'set dic.ole and of a 
current disc. 
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DATA r I IT INC 


Both the Non-linear and Spherical Harmonic Models V'^re fit to the 
data by the means of a least squares technique. This consisted of mini- 


mizing in both cases. Define as 


= y y 


j=li=l °i' 


where c. is the uncertainty in the i— set of measurements, , 62^-, 
are the respective x, y, and z components of the model B-field, and 
^li’ ^2i’ ^3i respective x, y, and z components of the experimental 

B-field. 

The expected value of N data points and with n parameters in the 

model is 3 N-n. In our analyses we accepted this value as being essentially 
correct, and adjusted a- so that a correct value of x^ was obtained. 

Uncertainty - We assumed the variance in each measurement was equal to 
the sum of a known variance due to digitizing error (see Instrumentation 
Uncertainties) and a second source of error due to uncertainty in the space- 
craft's position, instrument noise, time variations in the field, and 
probably other sources also. Although this second contribution to the error 
could not be determined experimentally, expected it to be roughly pro- 
portional to the magnitude of B at each point. Thus 

ct - . a B. 

where is the standard deviation in measurements due to digitizing 
error. The parameter a was adjusted until an aoprcpriate value for 


was obtained. 



/ii mrm zation rccjtmes 


Two different types of least-squares fitting were used in the two differen 
models. The first method is exact, but only works if B is a linear combina- 
tion of the coefficients which are being fitted. The second method is an 
iterative method which vie used to fit the dual dipole model, where the 
field was not a linear combination of the parameters vare fitting. 

Li near - If is a linear combination of the parameters to be fit, 
minimizing >;^ is fairly straightforward. (See Mathews & Walker, pp. 391-2). 

Let n 

E . . = 7 C . . A 

Ji Jim m 

where a^ are the parameters to be fit. Define a data vector, X, and a 
measurement matrix M as follov/s 


X = 


3 N 

I > 

j=l i=l 


Tim 


B'.-. 

Ji 


3 N 


M = M 
mt 


/.m 


C.. C... 

,nm Tlx 


a . ■ 
1 


j=l i=l 

The vector of the parameters, A is given by 
A = M“iX 

The value c^, the variance in the oarameter is oiven bv 

n ' ' m 


Quadratic- If B.... is not a linear combination of the oarameters, a_, 
j 1 ’ m ’ 

a more sophisticated approach must be used to find the minimum of We 

used a routine which provides second order convergence which was developed by 
Davidon (1966) and modified by Decker (1976). 

Goodness of Fit- The nocdness of each fit was found bv ccmoutinc the 


oercen 


r DMC 


(root mean souare) deviaticn oetween the model and sxoerimental 


fields . 




It should be noted in conclusion that the interpretation of least 
squares curve fitting needs to be done very carefully. In the cases of 
both our models the fits we obtained are sensitive to which data points 
are included in the trajectory, and which are not. In the case of Pioneer 
10, for instance, the trajectory included a latitude range of only ± 13°. 
In the higher order spherical harmonic fits, a very good fit was possible 
which included a large quadrupole and octupole terms. If the resultant 
fit is compared to the field along the Pioneer 11 trajectory, there is a 
good agreement with that part of the model which is close to the Pioneer 
10 trajectory, but a very poor fit in other regions. 
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DATA 


The data used for these runs was 1-minute averages from Pioneer 10 
and 11. The Pioneer 10 data was obtained between 2.8 and 6.5 Rj. It 
included a latitude range of ± 13° and a longitude range of 33° to 206°. 

The Pioneer 11 data was obtained from 1.3 R, to 7.0 R,, 32° to -50° 
latitude and 0° to 360° longitude. Unfortunately parts of both trajectorie 
included data which appeared to have anomolies in it. Both of these 
anamolous segments occured during occultation by Jupiter and resulted from 
an uncertainty in the spacecraft orientation. Unless noted otherwise, 
data was fitted with this occultation data removed. 



RESULTS 


The following results were obtained for the non-linear model: 


Model : 


Pioneer 10 Pioneer 11 (R>2.SRj) 


MAIN DIPOLE 


M (Gauss -R,“) 

3.709 

- 

.008 

3.902 r .005 

Lat, ^ 

80.71 


.04 

SO. 79 ± .07 

Lon. 

120.3 

r 

.6 

145.6 z .7 

Offset (Ri) 
C ^ 

-0.105 


.002 

-.11^ z .002 

C 

0.009 


.002 

-.004 z .001 

'z 

0.095 

t 

.003 

-.018 z .001 


RING 


Radius (Rj) 

7.4 

t .04 

8.0 

z .13 

Current (x 10° amps) 

.021 

z .001 

.225 

z .003 

Lat. 

81.0 

z .2 

68. 

z 1. 

Lon. 

217. 

± 3. 

125. 

z 4. 


SECOND DIPOLE 


M (Gauss -R,3) 

.292 z .006 

.306 z .006 

Lat. 

65.6 z .3 

35.8 z .5 

Lon. 

164.5 z 2.2 

86. ± 1.4 

Offset (R.) 



C 

-.75 z .02 

-.61 z .03 


-.27 z .01 

-.44 z .02 


-.82 z .02 

.29 z .01 


% RMS Deviation .449 .49'! 


For comparison, Pioneer 10 data v/ith just one dipole and a ring, and 


with just one dipole are reported below 



6.0 



Pioneer 10 x^i " Single Dipole 




'2 1 


* 7 “ 



Inholind 


Outbound 


Radius (Rj) 


[■ inure 1 




7.0 


f).0 


5.0 

InOound 


inv^ Dual Dipole witl> Rio9 
Pioneeir lO^^v ouai i ... 



4 « 6.0 6.0 7.0 

3.0 

Outbound 

Radius (Rj) 

Fifiure 7 . 



Dicole 


Rina 


M(Gauss -R,2) 

4.25 


.01 

3.936 


.009 

Lat. ^ 

79.22 


.04 

79.18 

- 

.09 

Lon. 

140.13 


.16 

134.0 

4 . 

.3 

Offset (Rj 
^x " 

-.110 

j. 

.002 

-.153 


.004 

Cv 

- .065 


.003 

.006 


.003 


.022 


.002 

.031 

r 

.004 


Radius (R,) 

8.2 = .1 

• • 

Current 10- amps) 

.079 r .006 

-- 

Lat. 

81.8 ± .4 

-- 

Lon. 

172.3 = 2. 

-- 


"c RMS Deviation 


1.085 


2.675 


Figure 1 below shows the values of x^^i ^ function of the radius of 
Pioneer 10. The large values of x^ near the end points indicate a systematic 
deficiency due to exclusion of the ring portion of the rrodel, since the relative 
effect of the ring would be largest at larger r. 

Figure 2 figure shows how this effect has been eliminated in the dual 
dipole fit which includes the ring. 


Spherical Harmonic Fit 

Different numbers of interior and exterior Doles were tried using the 
spherical harmonic approach to find the besi description of the magnetic field. 
The results for separate runs of Pioneer 10 and 11 data are summarized in 
Table I. 

.Notice that the fits to Pioneer 10 data with mors than two interior poles 
qivG values for q‘.‘ and whicn are obvious ly inconsistent with the other fiis 
as veil as with the non-linear model. We believe this is due ~c sne limited 
latitude range in the Pioneer 10 trajectory. 



SPHERICAL HARMONICS 





Pioneer 

10 


Pioneer 

11 



Pioneer 11 

R>2.8Rj 




lilt. 

Poles 

3 

2 

2 

3 

3 

3 

3 

2 

3 

2 

2 

3 

2 

3 

3 

Ext. 

Poles 

3 

1 

2 

1 

2 

1 

2 

3 

3 

1 

2 

3 

3 

2 

1 

% RMS Dev. 

.3356 

.8361 

.4763 

.4530 

.4006 

.7023 

.6330 

.9393 

.4500 

.9527 

.6541 

.3225 

.4201 

,3979 

.4823 


G10 

3.918 

4.286 

4.244 

3.057 

3.215 

4.080 

4.091 

4.025 

4.1274 

4.0277 

4.0283 

4.0989 

4.018 

4.144 

4.104 


GI1 

-.9082 

-.5712 

-.5522 

-.2511 

-.2983 

-.4505 

-.4398 

-.4887 

-.4325 

-.5069 

-.47808 

-.4252 

-.4509 . 

-.4653 

-.4697 


G20 

.3/91 

.2310 

.08604 

.02466 

-.1958 

-.1059 

-.1254 

.06419 

-.0920 

.11474 

.19540 

-.2132 

.1966 

-.07854 

-.04407 


G21 

-1.178 

-.7107 

-.5979 

-1.942 

-1.674 

-.8886 

-.8545 

-.9848 

-.726 

-1.0638 

-1.0519 

-.8125 - 

1.0204 

-.6932 

-.8148 


GI7 

-.2320 

.2464 

.3259 

.6183 

.6344 

.3463 

.3498 

.4034 

.3467 

.3247 

.41275 

.3142 

.5186 

.3865 

.3813 


G30 

.4996 



1.681 

1.193 

-.3133 

-.2992 


-.0535 



-.1730 


t.8089 

-.6205 


G.31 

-1.315 



.5093 

-.00495 

-.7620 

-.8303 


-.71792 



-1.0635 


-.6116 

-.5317 

Gduss -J 

GJ2 

.8220 



1.884 

2.011 

.3540 

.4105 


.3853 



.3735 


.5581 

. 3998 


G33 

.2054 



-.6966 

-.5417 

-.1228 

-.1310 


.38532 



-.2172 


-.03085 

-.04654 


1111 

.4258 

.5293 

.5470 

.5854 

.5739 

.6025 

.5996 

.4996 

.5793 

.61751 

.59228 

.5987 

.5196 

.5821 

.6018 


H2I 

.5173 

-.5760 

-.4720 

2.661 

2.356 

-.07034 ■ 

-.07621 

-.04309 

-.06945 

-.08030 

-.06207 

-.0623 - 

.04682 

-.007371 

-.03650 


1122 

.4458 

-.1406 

-.08553 

-.7880 

-.7194 

-.3534 

-.3142 

-.5810 

-.23897 

-.29898 

-.29405 

-.3049 

-.4263 

-.5073 

-.4300 


nil 

-.9805 



.06734 

.3483 

-.2811 

-.2587 


-.06567 



-.2770 


-.5155 

-.4725 


1132 

1 .007 



2.060 

1.909 

-.1840 

-.1870 


-.11097 



-.1702 


-.01961 

-.0807 


.1133 

.8oG9 



-.4213 

-.4733 

.2521 

.3102 


.29365 



.0892 


-.08441 

+.03338 


GIO 

65.63 

-181 .5 

-170.2 

189.0 

139.9 

-116.8 

-122.0 

-289.4 

-316.9 

-98.7 

-107.5 

-231.2 

-246.0 

-106 

-102.9 


GI1 

310.3 

25.22 

23.79 

-53.59 

-48.51 

12.65 

4.367 

43.56 

8.5 

29.8 

9.4 

-3.5 

-1.2 

11.34 

14.2 


G20 

-6.737 


7.184 


-3.821 


.9166 

-5.235 

5.2 


-4.0 

4.4 

-7.5 

.5982 



G2I 

4.184 


1.544 


2.089 


.9993 

-2.615 

-3.8 


1.2 

-.7 

-.3 

-1 .65 



G22 

.7551 


1.416 


-1.061 


.7723 

-1.742 

-.7 


1.6 

-.9 

-1.5 

.294 



G iO 

5.467 







1 .904 

-1.0 



- .3 

1.6 




G3I 

4.333 







-.4251 

-.2 



-.4 

-.1 



Y — 

G32 

-1.446 







.9232 

1.3 



.4 

.5 




G33 

.006711 







-.0737 

.4 



.4 

.2 




1111 

35.08 

-10.46 

-1.866 

-28.05 

-14.28 

-49'. 32 

-39.68 

39.04 

-38.7 

-63.2 

-61.1 

-47.5 

15.7 

-28.93 

-42.28 


ii.’l 

-6.730 


3.899 


-.5372 


-.1775 

-1.063 

-4.9 


2.4 

-4.0 

-.7 

.4459 



il22 

-6.926 


.1672 


.7495 


-2.744 

6.737 

.2 


-2.4 

-1.0 

2.0 

.2114 



iin 

.9/74 







1.205 

2.9 



1.4 

.9 




ll 12 

-.6201 







-.2174 

.2 



.5 

.2 




Iil3 

.1719 







-.4314 

-.8 



-.3 

-.1 




Table I 



Table II shows haw the values calculated in the soherical harmonic 
exoansions relate to the paraneters in the other model. 

The coefficients don't agree between the various fits in the table for a 
number of reasons. 

1. We believe the higher order fits of Pioneer 10 data are not good 
because of the small latitude range in the data. The best fits of Pioneer 10 
data are most likely the 2,1 and 2,2 fits. 

2. Adding additional terms to a truncated series will affect the least 
squares fit of all of the terms in the series. We would expect the best values 
for g? to be in the fits using three interior coefficients. 

3. There may be seme secular variations in Jupiter's magnetic fields v/hich 
caused a difference between the Pioneer 10 and Pioneer 11 coefficients. 

For the above three reasons we decided the best description of Jupiter's 
average field would be provided by a fit involving three interior and three 
exterior poles and using combined Pioneer 10 and Pioneer 11 data. We included 
the exterior octupole coefficients to lend the necessary significance to the 
exterior dipole coefficients. Because of the axial syrmetry of a disc field, 
we would expect the exterior quadrupole terms to be zero. Hence, it is 
necessary to include the octupole terms to avoid generating an asymmetry in the 
series by the way we've truncated it. The results of this fit are in the 
following table. 
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SPHERICAL HARMONICS 


OJ 

c:j 



M 

Gauss 

4.356 

4.314 

3.280 

3.123 

4.149 

4.174 

4.210 

4.158 

4.085 

4.077 

4.100 

4.106 

4.190 

4.044 

4.164 


Dipole 


Eat 

Lon 

CX 

"j 

CY 

79.71 

~1567r~ 

-.096 

o 

CO 

o 

1 

79.6 

135.3 

CO 

o 

1 

-.068 

78.6 

117.5 

-.319 

.395 

78.2 

113.2 

-.383 

.460 

79.6 

126.8 

-.135 

-.010 

79.5 

128.0 

-.125 

-.006 

79.8 

128.6 

-.109 

-.000 

79.7 

126.3 

-.129 

-.011 

80.1 

134.4 

-.152 

-.005 

80.3 

131 .0 

-.155 

-.013 

79.3 

128.9 

-.155 

-.017 

78.8 

129.4 

-.156 

-.016 

80.1 

126.7 

-.109 

-.012 

75.6 

154.9 

-.134 

.054 

79.8 

125.4 

-.123 

-.008 



The computation of J assumes ana of 1.6. 

In/OiitiB^Both I “In O^Out 
()cc?:N-No Y-Yes C-Close data exluded 


Table II 



Interior Coefficients 
(Gauss 1 

= 4.C99 ± .004 

** 1 

gi = -.486 ± .003 

gO = .034 r .008 

g* = -.81 ; .01 

g^ = .303 ± . 007 

g_” = -.04 ± .02 

gi = -.22 = .03 

g= = .41 ± .02 

3 

g- = -.45 = .02 

= .526 ± .003 

1 

ni = -.084 : .007 
2 

h2 = -.404 ± .006 
2 

hi = -.04 ± .02 

3 

= .22 - .02 

hS = 0.02 = .02 


Lx-erior Coefficients 
f ) 

; i 

0° = -121 ± 5 

gl = 9 r 2 

1 

g° = - 1.1 r .2 
9: = -.4 ± .2 

= -.1 - .2 
g= = .13 r .0^ 

gi = -.35 r .04 

3 

a- = .12 = .05 

" 3 

g- = -.05 = .03 

Ri = -17 ± 4 

1 

Ri = 2.1 = .2 

F 2 = .1 = .2 

2 

B1 = -.30 r .07 
Ft = .27 : .04 

fi' = -.24 ± .03 


These correspond to an offset dipole of r^agnitude 4.16 Gauss-Rj^ with its 
axis at latitude 80.1' = .90, longitude 132.7° z .4°. It corresponds to a ring 
with = S3 y, with its axis at 80.7° t 2.3° latitude, 118 ± 10° loncitude. 

The percent RMS deviation was 1.1 %. With = S3 y, the current density 
would be .06 amps/meter. Assuming the ring goes from 10 Rj to SO Rj this 
would mean a total current of 7 x 10^ amps. 
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The wo different approaches show excellent anreefnent in the magnitude of 
the total vector dipole moment. The two values obtained with the nonlinear 
model vere 3.94 and 4.166 Gauss-R,^. The value obtained in the soherical harmonic 

U 

approach was 4.16 Gauss-R,-. The dipole is tilted 9.3°, 9.2° and 9.7° from 

U 

the spin axis in the three cases. The .longitude of the dipole is 126.3°, 145.6°, 
and 132.7°. The larger difference in the longitude values is due to the high 
latitude of the dipole axis, as v/ell as the lack of high latitude data on either 
mission, especially Pioneer 10. 

The placement of the second dipole or the exact nature of the external 
field sources is more uncertain. The disc current tends to roughly line up 
with the main dipole field. The fit given by the spherical harmonic analysis 
agrees fairly well with the fit obtained using just a single dipole and a 
ring with Pioneer 10 data. This goes well with our intuitive feeling that the 
disc ought to line up fairly well vnth the main dipole 'ield v/hicn is containing 
the current. 

The fact that we vere unable to get a consistent fit for the second dipole 
is probably an indication that the interior field is more ccmplex than was 
assumed in our simple dual -dipole aporoach. Further confirmation of this is 
the fact that we v/ere unable to find a good fit to the data if Pioneer 11 data 
with R < 2.8 was included in the runs. The location of the second dicoie 

U 

may be further complicated by secular variations in the field sources within 
Jupiter. It should be noted in passing, however, that the locacion of the 


second dipole in both fits is roughly at the same longitude as the field 


ancmolies reported by Serge and Gulkis (1975), Conway and Stannard (1972,1975; 



and ctners. 



An analysis of systematic errors in the fits of the soherical harmonic anc 
dual dipole models reveal two interesting things about tne model. The first 
conclusion is that the current disc probably starts within 9 Rj of Jupiter. 

There are significant increases in the percent field deviation from the linear 
model between 6 and 7 Rj as would be expected in a truncated spherical harmonic 
expansion of the field. It would take many additional higher order terms to 
generate the larger increases in B one would expect near the edge of a current 
disc. 

A second systematic anomoly in the data was noticed in both the dual 
dipole and spherical hannonic fits. There was a bad fit to the data at about 
4.8 R, that was oarticularly noticeable on both Pioneer 10 inbound and outbound. 

U 

The effect was also seen, though not quite as pronounced, in the Pioneer 11 data. 
Such a bulge could be tne result of seme current flowing at about that radius 
since the initial assumptions in both the dual dipole and the spherical harmonic 
analyses are that there is no current within this region. If there is a 
current in this region, however, it would have to be a fairly small current 
compared to that in the outer ring. The current would also have to be such 
that the trajectory of both spacecrafts did not pass through it, since there 
is no evidence to support such a crossing in this region. 



V 


The oest average description of Jupiter's magnetic field is given by the 
data in Table . Physically, the field corresponds roughly to a field due to a 
dipole of moment M = 4.1 Gauss-RJ^ tilted 9.5° with respect to the spin axis 
and at longitude 135°. The external field source is a current disc with an 
inner radius between 7 and 9 Rj and with an axis that roughly corresponds to 
the axis of the main dipole. 

This type of analysis v/ill permit better resolution of tne various fiel 
sources as more field data is gathered during future Jupiter missions. 
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riONLIilK 10 MUASIJKI£MHNTS OF 
JUm i-R’S MAONFTIC FIFLI) 

I), li. Jones,' H. J. Sinilli,^ L. Davis, Jr.,* 
I). ,S. ('ollnirii,'' I*. J. ('oleinaii, Jr.,* 

I’. Dyal,'* and ('. I*. SonneU* 


Alin u IliiJii of 642 days, I’ioiicei 10 passed the planci Jiipiler al a 
disiaiicc ol 2 K4 Jiipiler radii (KJ) al U2.U U.1'. (Local Jupiler) on 
hcccinliei 4. I‘)7.l. Diiiiiig llie period Novciiiliet .JUlhiougli Decemher I2 
a vector lieliiim inagnelomcicr oblaiiied iiieasiiicmcnls of llic Jovian 
iiiagiielic field The purpose of lliis paper is lo provide llie pieliiiiinary 
lesiilis of lliese iiieasiircineiils. 

I’loiieer ID appioaclicd Jupiler al alioiil 6.5° soiilli plaiieloccnliic lali- 
liide (Noveiulicr 30) al a .Suii-Jupilcr-Pioiicer angle of 35°. A single how 
shock crossing was ohserved inhouiul al I OH KJ (Figure I). The inagnelic 


lO*. 


lO'* - 


PIONEER lO ENCOUNTER 

JOVIAN FIELD MAGNITUDE 

. _ j ^ j j y I |- 

PCRIAPSiS 

9 

t 
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I'lgiiic I: I'lve-iiiiiiiile averages of the held iiiugiiiliule llirougli periapsis. 


I llrigliaiii Young lliliversily, Provo, lU' 

2 . il l I’lopnision I alioialoiy, Pasadena, ('A 

3 . ( jIiIoiiii.i liisliliile of I'ecliiiology, Pasadena, CA 

4 . Ames Uesearch Ceiller, Mollell Field. CA 

5. llniveisily ol ('.ilifornia, I os Angeles. CA 
It. Univeisily ol Ail/on.i, liicson. A/. 


field jumped from 0.5 to l.5y (10'* gauss), and the magneioslicaih fields 
hehind llie shock vaiied iiiegularly in magnilude and direciion. Waves were 
ohserved in Ihe field piupagaling upstream prior lo the crossing of llie 
shock. Ihe inlerplanclary field direciion oiilsidc the shock was such ihal 
energetic pailicles could propagate upstream lo Ihe spacecraft as ohserved 
hy Ihe radiation deleclors. 

The niagnelopaiisc was ohserved at 96 KJ, which was earlier lhan 
expected based upon a simple scaling of Ihe corresponding l-iailh geoineiry 
(.see Figure 2). Ihis implies either a standoff distance lhal is lelalivciy 
small as compared to l-aiih oi else an outward motion of the magneto- 
sphere. Al Ihe magnciopausc Ihe Field jumped ahruplly to Sy. Ihe cor- 
responding magnetic energy density just inside the magneiospherc (I O '" 
eigs/cm ') would appear lo he insufFicienl lo withstand the pressure of a 
nominal shocked solar wind (estimated lo be 5 x 10'"* crgs/cm* ), or cl.se 
the plasma density was much less than the nominal value of 0.2 cm *. If 
the solar wind density was nominal, then ihis implies lhal ihe piincipal 
magnetosphere pressure was due to a fi^4 plasma inside ihe magnelo- 
spliere. A lower solar wind density implies an outward tnolion of the 
magnelosphere. 

J'he Field inside the magnelosphere exhibited a persistent southward 
componenl. Hence, the field lines were probably closed and ihe orienla- 
lion of till; dipole source at the planet was roughly parallel to Jupiter's 
spin axis, as inferred from radio astronomy measuremenis. The Field 
magnitude remained near 5 t from 90 KJ to about 50 KJ, hut was very 
irregular, with freqiieni dips to It or below occurring. 

riie Field in the outer magnetosphere was strongly distended such lhal 
ils direciion was elongated parallel lo Ihe equator. Uiere was no well- 
defined orientation of the Field into magnetic meridian planes. Keferring 
the field vector lo a Solar-Jiipiler (SJ) coordinate system (X = S, S from 
Jupiler towards the sun; Y = normali/ed J X S, J parallel to Ihe spin axis 
of Jupiler; and Z compicles ihe right-handed system) il was noted that 
much of Ihe lime the X and Y components were of opposite polarity, 
suggesting a spiraling of the meridional planes of I) due lo plasma effects 
causing the fields in ihe ouler magnetosphere to lag behind those closer to 
the rapidly loialing planet. There was a frerpieni interchange of polarities, 
and occasionally there were periods when they were the same polarity. 
The X componenl was usually negative and the Y componenl usually 
positive. Hence, the spacecraft was prohahly below the symmetry plane 
most of the lime. The occasional coupled X and Y polarity reversals sug- 
gest that al limes the symmetry plane passed underneath the spacecraft. 
These results and Ihe ahnormally thin inbound magncloshealh an; con- 
sistent with the existence of a relatively Hal, “disclike” outer magiiclo- 
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I'lnuic 1 . I rii)ccii)fy |)lol of the Siiii-Ju|iilcr-l'ioiu.'cr angle stiowing bow 
sli()i.'k and iiiagiieto|iuii!.c tro!, sings. Also shown are llie rclalive posilions of 
the ('uiii:iil slii'ol ciossiiigs. 


splicie. Sucli .1 ilisclikc field requires llie presence of a ciirieiit slicel 
miiglily in ilic syniincliy pl.ine, and il is teinpling to infer llial llie occa- 
sional dips m llic field inagiiiMide to < I 7 may have resiilled from llie up 
and dow'ii movcincnl of ibis slieel curreiil pa.sl rioncer 10 , peiliaps in 
icsponsc lo t lianges in llie solar wind, cic. 

flic licid slicngtb began lo lise monoionically al aboni 25 KJ, and 
pciiodic vaiialions in field diieclioii baving a lO-lionr period were luiled, 
indicating cnliy inlo ibe inner magnelospliere. Sliorlly ibeieailer, periodic 
cllccls in field magnitude became discernible, wbicli appeared lo correlate 
well wiib ibe ebanging niagnelic taliinde of 1 ‘ioneer based upon ibe 
nominal ladio asirononiy values for the longiinde and iiiclinalion of the 
1101 lb iiiagnclic pole. The maximnni field slreiiglb measnied was O.IK 
g.niss, coiiespoiidmg lo a niagnelic nuimeni of aboni 4 ganss-RJ^. 

I be passage oiiiboniid (brongli llie magnelospliere was al a Snn- 
Jnpiici-rioiicci .ingic ol lOO’ al 05.10 bonrs local lime (bigiire 2). Ibe 
InbI was mneh inoie icgnlar lliaii dining llie inbound porlion, with clear 
evidence ol Kl boiii peiiodicilies onl lo *)() KJ. As llie lailial dislaiice 
iinie.ised lowaids ‘>(1 KJ, llie field again became slioiigly extended 
piincip.dly III .1 ladial iliieclion, so llial il lended to lie paiallel lo ibe 
eqii.iioii.il pl.nie. Ibe field lay nearly in Ibe local meiidian pl.me in Ibe 
innei in.igneluspbeie, bnl a syslemalic devialion of lens <d' degrees 



BYKJ 



Figure 3; A portion of tbc niagnelic field data of day 342. Tbe observed 
ebanges in magniliide anil direction are consistent with Ibe passage ol Ibe 
spacecraft inlo a current sbeel. 
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rii>iiic •!: i‘luiicloci:iilrik: U-I.iititiule plul of llie cstiiiulcil poNilion uiul 
extent til e.n:li eiiiieiit slieci Cfiixsing. 


(lcvelti|U'il liiwnuh llie aiilisiilar iliiccliim, a|;aiii coii:ihli:iil will) tlic 
app.iieiil xpiialliiig of llic fieUI <icen inhouiul. 

Dlpx ill llie licltl slreiigtli willi occasional partial reversals in the ndd 
ciiniponents ot;cinreJ at lO-lioiir intervals similar lo the one displayed in 
I igine I (see also l-'igiiie 2). The plancloccniric (roialing) cooidinalc 
posiiion of the esiinialed centroids and laliuidinal and longiliidinal exienis 
of iliese tlips aie ilisplayed in l-igiires 4 ami 5. l-ignie 6 shows an idealized 
lepieseniaiion ol ihc enrreni sheet required lo explain ihe field dips and 
leveisals nhseiveil ihiiing the oiiihoiind passage. 1‘ioneer cio.s.seil ihe inag- 
nclopaiise on the oiilhoinul pass at 98 ItJ, a localion consisleni with a 
cylimliii.ally syniinelric, disclike shape for ihe niagnelosphere near the 
etpialoiial lei'inn ( Hgiiie 2). 

OhseivaiiuMs of Jiipiiei's magnelic field in Ihe radial range 2.84 lofi.O 
KJ wcie used lo ohiain a hesi least stpiares fil planetary dipole. During (his 
poilion ol ihe lra|ecloiy, llie standard plaiielocenliic (J(i) laliinde 
(I.A I Kit ami longiimle (I.ONJG) varied from I.V Ui i|.t“ and 179" C'W 
lo do" (ol I8l“ lo .114“ in Syslem III) respectively. The coires|>onding 
ilipoir IS I haiai leii/ed as I'ollows: 

Miiiiiciil 

M = 1.867 Gans^ KJj 
Ml ATJG ^78.69° 

Ml ONJG ^ I1.S.18" 

= 224.42“ Syslem III 



I'igutc 5; IManeloccnliic K-I.oiigilnde plot of Ihe estimated centroid 
position of each current slieel crossing and the estimated average variation 
of the longitudinal extent of the tiaversal into Ihe sheet. 
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I il^iiii- o: All iiK-.ili/cil syminctiic current slieel imnlel int'errcil Iroin lliu 
>>Iim:ivc>I oiilliuiiml licKI d.ila. 


Offul 

( X ^ -.N5 RJ 
( Y I O.R) RJ 
(7. ^ 1.070 RJ 

VARMXVZ - I .OH X 10 ) (laiiss^ 

VARItXY/. - I 6.1 X 10* CaiissJ 

wlicic VARM.XS / is ilie iins of ilie compoiietil vaiiaiices of M, cic. (Noie 
ili.ii iliese v.ihicN ilillci sliyjilly I'liini lliosc re|iiiileil eailicr by .Smilli, el al., 
I‘»M) 

III llic i.iii)’c limii J.K4 Id 0.0 RJ, llic i>ll'ceiller. Oiled dipiile derived 
.iliDve .Mill liy .Siiiiili, cl al. (l')7-t), Ills llie Pimicer 10 inagnelomelcr dala 
■ iliscivcd Dvci Oils laiif'.e In williiil alioiil of llie Held, file iiiaxiniiiiil 
siiil.u'c field pu'dicicd by ibis model is nearly 12 Caiiss, which is coii- 
sisiciii willi llic iiiaxiiiiiiiii Odd sireiigth leipiired by ihe ohseived dcca- 
lilclilC l.idl.llioii. 
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A STUDY OF SINGLU AND DUAL DIPOLF 
MAGNETIC FII-LD MODELS FOR 
JUIMTEU: PIONEER 10 

Douglas IL Jones and John G. Melville 

Hrigliani Young University 

Hie lestills of a iicutly real lime analysis of Ihe magnetic field data from 
Piiiiieei 10 have hecn reported previously (Smith, el al., 1974). 1'hc 
|)iii|»ise of this note is to report an analysis of the near-Jupiter data 
covciiiig roughly the same range of trajectory parameters hut using as 
much of the data in this interval as possible. In addition to recomputing a 
licsi t paiamcler fri (rmninmm variance of M), we also computed best 
(i parameter and 1 2-parameter fils (minimum variance of II) to the tiala. In 
the I 2 paiamcicr, or dual, dipole model, we have followed a suggestion by 
Conway and Slaimard (1 972) which they proposed to explain an anomaly 
seen in ibe radio astronomy linear ptdari/alion data al approximately 
System III longilirde 22U . Tbe second dipole is postulated to lie near lire 
sill lace ol the planet. 

'Ibe two basic tecbniipies used to obtain the best model parameters are 
related to Ibe well-known equation for the magnetic field of an offset 
dipole. 



wlieie KC - R - C 

It - position vector of Pioneer relative to a planet fixed 
(loialiiig) coordinate system 
C - ollsel of dipide 
I his can he wiilien in maliix form as 

(ll)MA)(M) (2) 

Ol . .ilteiiialively 

(M) = (A)' (II) (3) 

whi-ie 


Uil 
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(RCX)^ - 2/.l(ItC)^ (RCX)(RCY) (RCX)(KCZ) 

(RCX)(RCY) (RCY)» - 2/3(RC)2 (RCY) (RCZ) < 

(RCX) (RCZ) (RCY) (RCZ) (RCZ)J ■ 2/3(K(’)^ 

In Ibe M method, the dipole offset C is varied (i.e., 3 parameters) until the 
rrns M variance is minimi/ed using equation (3) above: The II method 
requires that both M and C be varied (i.e., 6 parameters) until ibe rrns II 
variance is minimized using erpialion (2). If the measured field is due to a 
simple olfsci dipole, tbe two methods should give nearly the same “best 
ill" dipole. However, if the disagreement is outside of the ex|)ccled experi- 
mental error, then the source cannot be a simple dipole and tbe (i- 
paramcler model should be the best. The rrns R variance for each model 
can be used as a quanlilalivc means of determining the best model. 

In computing the variances, we ftrsl transformed lire measured II from a 
Pioneer Inertial system (Pli) into a planet fixed rotating coordinate system 
(JG). Interpolation of the trajectory data accurate to at least four decimal 
places was used to obtain Ibe values of R corresponding to the midlimc of 
each data sample. Tbe observations of Jupiter’s magnetic Held over tbe 
approximate radial range 2. 84 to 6.0 Jovian radii (RJ), the System III 
latitude range -I3" to -t-IJ", and longitude range IHF to 3I4° were used 
in tbe present analysis as before. 

Using the M method, the single offset dipole that gives the best least 
srpiares 3-patameler fit to the Pioneer lO magnetometer data used has the 
following parameters: 

M = 3.867 Gauss RJ3 
MI.AUG = 78.69“ 

MI.ONJG= I.35.I8“ 
ex = -.145 RJ 
CY = + .030 RJ 
CZ = I .070 RJ 
rrns var M = 1.08 X 10 ^ 
rrns var II = 4.(>.l X 10 *’ 

The ntajtrr ditference between ibis dipole and that reported earlier 
(Sirtilli, cl. al., 1974) re.sulls frtuu tbe use of ten-minute averages which 
have been collected for an ciior in the roll alliiiide of the spacecraft and 
fur an electronic phase lag occurring in the iiiagnelumeler. 

1'lic variable metric minimi/aiion program developed by Davidon 
(1966) was used in the model studies based upon the inis variance of II 
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(I e., ilic 6- ami 1 2-paramcler iiuulcis). The lies! leasl s(|iiaics fil single 
iill^L-i (li|iiile iiuulcl (ihlaincd using this iiielhod has ihe paiaineleis: 

M=4.07 Gauss 
MI.AU(; = 78.61*' 

MK)NJG= 139.8*“ 
ex = -.125 UJ 
( Y = -.038 RJ 
(•/. = + .05 1 RJ 
iinsvar 11= 1.05 X 10 ‘ 

r'lilliiwiiig ihe suggeslion hy Conway and Slannard (1972), we have 
im it'.i>.cd ihe degiee of eoniplexiiy ol Ihe models hy allowing lor a second 
dipole. Ihe dual dipole model lhal Ills the data hesi has Ihe paiameters: 

l-'inl Ihpiilf 

M =4.022 Gauss RJJ 
Ml AUG =80.08*“ 

Ml ONIG = I 33.(.0'“ 

CX = -.200 R j 
CY = - .03386 RJ 
(•/. = + .06812 RJ 

SiCKiiil Ihpiilc 

M = 0655 Gauss RJ’ 

Ml. AUG = 12.5.1*“ 

MIONJC. = 1 41. .5*“ 

CX = - .8221 RJ 
CY = i .1600 RJ 
CZ = - .1184 RJ 
imsvai ll = 7.')7 X 10 '’ 

Om- miles llui ihe second dipole is al ahoiil .95 RJ, and (he offset in the 
eipialoii.d plane isal ahoul Syslem III longilude 20*).3'“. 

Iloili of ihc II models fil ihe data siguificanlly heller lhan does ihe first 
oi M model. Ilowevei, the model which gives Ihe lies! Hi iscleaily Ihc one 
whii h im hides a secotui dipole. One notes lhal Ihe tills of all of Ihe 
models agiee leasonahly well with values derived fiom ohservalions of the 
vaii.ilioii in ladio inlensily (llerge I'l7.t) and liom ohservalions of 
Jupiiei's decimeiiic ladialion (Roheils and KomesaioH l‘i(,S ). | he dual 
tlipide model also appears lo he consistent with an inleipielalioii of Ihc 
pol.ui/aiion anom.dy suggested hy Conway and .Siann.ud (|U72). Uiis 
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model also has a Syslem III longilude for Ihe main dipole which is mosi 
consistent with predictions based upon tlie radio mcasuremenis (llerge 
1973), although there is still an apparent discrepancy of about 3*“, or 
1.3%. 

It is difficult lo justify the piesence of a dipole source so near the 
surface of the planet. Other more complicated, but physically more 
reasonable, field configurations for the anomaly are possible, but most aie 
far mure difficult to work with in studies of this kind. If subseipieni 
observations by Pioneer 1 1 confirm the existence of an interior field 
source near Ihc surface of ihe planet, considerable information regarding 
the dynamo origin of planetary fields in general, and the struclure of ihe 
interior of Jupiter in particular, will be obtained from these data. 
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from ll»c rnirfi\Iiiig% of the Vlah Academy of Scieiicex, Arts, and 
Letters, vol. 52, pari I, 1075 

I’KI- UMINAKY MODIIL S l UmiiS 01- Will 

ma(;ni:k)si'iikrlof jumi-u: I'lONi-i-it lo 

Doiiplas I:. Junes :iiul John (■. Melville 
liiiglmm I'oinig Ihiieersity 

liilrotliiclioii 

'I'lic I’ioiieer 10 spacecnift pii.sscd witliiii 2.8*1 Jovian radii (l\j) of (he 
|)iaiu't jiipilcr on Dcer-inlter *1, 107.1. Extensive ohservations of the Jovian 
iii.i)pit'(ie held and i(.s inleraelion with (he .solar wind |)laMiia were made 
while (h<! spaeeerali was within alMiiit l(K) ily rrf the planet, 'lire iiiaf'neto- 
sphere was loinid to he S4*ver«*ly stretehed l>r*eanse ol the pres«*nee of 
an intense; enncnl sheet, which was partienlarly <-vident during the oiit- 
honiKl p.issai;e ol I’ioneer lO near the dawn ((‘iininator (.Smith, i*t al., 
I07’t). riots ol the an)'l<‘ between the orientation of the onthonnd held 
and the raihns veeloi hoin the planet to the S|)acecrart showed a strong 
teinleney lor the held lo become radial at large distaiu'cs from the plaiu't 
(sit |■'igln'e 8, .Smith, et al., 107-I). A similar trend has alsti been seen in 
both the inbound and ontbonnd IMoiieer II data (Smith, et al., I07.'i; 
Join's, el al., 1 07.')). \\’«* lepoil hei»‘ some pieliminaiy work on a nialhe- 
inalical model ol the magnelosphei'e of Jupiter whii'h is based upon the 
I'loneei II) onibonnd il.it. I. A pieliminaiy model study related to (he 
ontbonnd I’ioneer It) data lias al.so been reported by (ioeii/, el al.. 
(lOT'l). However, we have noted some Iniidameiilal conceptii.il errors 
in their sliidv, and it is al.so the [impose ol this paper to re|)or( a corree- 
lion III this e. II lie! analysis. We will also discuss some ol the implications 
ol the r.iili.il lield eonligiiration ililcrred I'roiii the l‘ioneei lO and 1 1 data. 

I'he .McIIkkI 

Smee it is .ilu ,ivs true that 

V ■ /I = 0. 

one can ri'|iiesenl It by an expression ol the liiriii 

/J - V/X Vg, 

where / .mil ;; are scalar Innelions of the eoordinales that are some- 
Imies leleiieil In as |■'lller potentials (l''nler, I7(i!); Trnesdell. l!).')l; .Stern. 
|!Hili). I'he iitililv ol this manner of repiesenting /i lies in the fact that 
Mill e yi is l.iiij',eiil lo the iiiterseelion ol the siiriaees / - i nnslaiit and 
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g =* consliiiit, (his allords a ilirect melliod of (.‘valiialing the shape of the 
lines of force. For exani[)le, for a tlipole in S[iherical cooidinati s we have 

f Afsin^O 

/““-p— ■ 

where/) = /V/i/>and 

«A* 

(In cylindrical coordinates, p will represent the dimensionless component 
of i( that is perjtendicnlar lo the axis of the di|M>le.) 'Ilie / function 
can be easily manipulated into the well-known constant L representation 
of a dipole llchl line, namely 

p—L sin^ 0 . 

Although the law of snpeqmsition holds for magnetic llehls, this is 
not generally true for the fnnctioiis / and g, i.e., 

VFx vr. = V /, ) X v( 

' 1 ' ' t ' I 

Alternatively, one notes (hat 

V/X Vg = V X (/Vg), 

.so that the vector potential A is related to the /anil g fnni lions (hrongh 

A =/Vg. 

For axi-synirnclric fields, / i.s independent of ^ and g = <l>. The vector 
potential in .spherical coordinates is then 


'■•liliif. a '*'• 

p sin 0 

and in cylindrical coonlinates, 

X,„ . .Ohlly 


For axi-syinnietric fields, or f functions sharing (be same g Innclion. 
one writes 


A= S/Vg, 

I 

or, allemalively. 

vx(s /Vg )- v(s/. )x 
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\Vc write willi Hi, heiiig the dipole fielil, 

a = il„ + il, + Hj. 

or 

Il - V(/„ hfi)X V«, + V/, X Vfi,, 

•SO (litil the |)ei ttiilKilion field ii,, iii given by 

lip => ll — il„ 

\ = il, + i/j 

= V/, X Vg, + Vf, X Vg,. 

ll I represeiit.s tlie axi-!>yiiiiM(!tric' portion of tlie pertiirhution field and 
if 2 the eoinpoiient eoiitrihiiting to the spiiuliiig. 


Spliericiil i’oliur Coordinate Model 


In Nplieiieal coordinates we have 

ji = ._L__J/lp_ 
* p^siiifl do 

ll2 - 

p sin 0 do 


Since tliere is clearly spiialing of the fiehl (Sinitli, et al., 1974), we have 
wi illen 


= «i + 

= + kp, 

where kp icpresents the spiraling, necaiise Hi and H,, share the same g 
Innction, ami therelore superposition holds for the respective f liiiu'' 

. tiinis, we will concenliale on thesi* coinponents of the iiehl only. 'I1ie i 

Innction /' -- f, L f], will then represent meridional plane projections 
ol the incasni c(l field. 

In di'iiving an fi fnni'tion, we start with a component ol thi‘ peiTiirha- 
lion Iiehl whose Innclional form may he easily deduced from the ilata. 

.Since the r.nlial eoinpoiient ol the field decreased and at times reversed, 
u'hich is I'onsisicnl ss'ilh passage into ii thin (.'iirrent sheet (Smith, et al., 

|!I7 I), a innclional foriii liir />, that is consistent with these lactors (see 
al.so Ihi (I. Ifir.'i) is 


I 
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C = 770, 


cos = 0.025. 


Alllioii^'li filling llic Fioiicer 10 oiiIIioiiikI data <|uiU‘ well, lliis F riiiie- 
lioii e.sliiliiled latlier aiionialuiis lieliavior at lilt'll laliltules. Since (lie 
(r.spressjoii lor />, can iiu:liide additional i'linctions of 0 wliicli are small 
near 0— n/2 (i.e., i'nnetions of cus 0). ono could write 

h = — tanli + c(is 0 li(p, 0). 

^ p“ costfo 

Kollowin^ lliis lead, alternate riinctions can be derived. One of several 
wliii'li fit the data reasonably well is 


/i= - 


A cos r cos 0 , 

-7.-2- [logr«sl»--; + 


cos 0^t 


, 7 ' J • 


wliere llie coiislaiits are tin* same as tliose listed above. 'i1ie correspond- 
ing; !■' riiiictjoii is plotted in l''i);nre 1 with Af = T X 10' (Smith, et al., 
1071; I!l7r>). Altlioii^b this function exhibits better behavior near the 
iiia^iietic axis, it is still iinsatisfaclory here. Heplaciiig Ce' ’’"'/p*' 
by — lo^; cosli(lA,-os I7„) produces an f, whicb matches the b,, data and is 
widl behaved at 0°, but iiisiiilicient .southward field results because 
ellects doe to iiia^iielopaiis<> clirrenis have not been included in the model. 
Clearly an additional f function is needed, but inlinite series technicpies 
will likely be reijiiired. 

Ne};le« tiii^ the presiiice of the magnetopause, the last closed field 
line crosses the iiia^iietic eipiator at p,. — MO for (he preceding models. 
Undi:r these conditions magnetic field lines originating at higher latitudes 
woiilil not cross ihe ixpiator and would thercliire be considered as being 
open. 

The Cylilulrieal Coordiiialv Model 


(ioiot/, et al. (Itl7-i) have ileveloped a model in cylindrical coordinates. 
l''oi this case the axi synmietric poilioii of the |)ertiirbation field is given 
by 



p dp 


MAGNETIC 

AXIS 



Ki^iiie I. Meiiiliiiii;il |il.iii<‘ ie|ireM-iiliiliiMI all llie iliii^iia-laas|alirrii' lii lal ,i\ alria'lii|iial 
ill a Nplii'iiaal aaiimliiiale ri'|iii'M'iilaliiiii. T1ie soliil anal al,e>lia'at aiaaiwA ii'|ira'Miil ilia- 
ava-ra^aa* lia'Ial aliia'a liiiii iiia'.asiira'al al sa-va-ral jiailiils alain^ ilia- aiaatlMiaaaial laa|a i laaiia". ail 
I’iaiiia a'i III aiail II, ri‘'|ia'i liva ly. Ilia* sliiaila'al liall aia^iilar w iallli |iaiilauai iia'.aa llaa* 
Ilia^lia-lia' a*a|iialaii ia'|iia*sa'iil\ .i |iaiiljaiii ait Ilia* a*a|lialaiiial ciiiia'iil ^lla a l a aiailli;iil.aliaiii 
asMiiliaal in lliis iiiikIi'I. anal a liisl airala-r a|i|iriiiiillialiam lai ilia- iii.ai;iia'la)|uaisa' 
iKiiiiialary is .ilsai iiialia-ala-al. I(a*)>iaiiis aivi*r wlifi-li llie liiiialiiiii Is la-li.ilili' aii- iiiilii ali-al 
III llie text. 'Ilia* ciiiva*s li*ava* Ilia* i< = 2 sjiliere til is|iially S|iaceal aii)>iilai iiileivals. 

wliere the clirreiil sheel hall-wiillh, P, could ba* soma* fiiiiailaan ail p. i.e., 

P = Piiplpu. 

for ii const, nil angular width sheet, or. in general. 


,‘\s beloie, a liinctional liiriii for b, (hat is consistent with (be ciirrenl sheet 
dala, etc., is 


P - Pii,if>lpi,)'‘- 

Since 



20 


in All ACADI MY I'lUK'l I niNl.S. VOI.. 52.rAKT I. 1175 


“ ~, 7- 1 

Tin: L'lirrcspomling liiiiclioii lor is tlicn 

Al)(,i - 1) hAz dC(p) 

-;r- '"B - ^771 

Note dial lor zID S 3, log cosh z/D is very nearly zJD — In 2, so that 




z Af)(« - 1) 


-■ In 2 + 


I’loUiiig vi'iMis p lor all zID^ 3, and h, versos p for fixed values of 
zin, allows one to determine the constants in (•oertz, et ul. (1U74) 
have plotted /;, in this manner and find 

«!■/<+ I = 2.77 


so that 


> 15 AD 

-zr 


/, = — — -r (log cosli zID f- 15). 


However, Ooeilz, et al. (l‘J74) have plotted b, versos p. where 

b. ° b^y/ 1 + bjb^ 


/;* = - kpb^. 

As a I (-suit, they ohtain a power law representation of the component 
paiallel to the magnetic (‘ijimtor which lies in the curved surface repre- 
sented hy 

0 -f Ap “ constant. 

liowevi i, the i( adtiiig/ will he for such surfaces, dm-s not represent an 
asi-symmeli le iield, and therefore cannot appropriately he added to the 
dipole J fimelion, which is axi-symmetric. That is, the / functions to 
he .uided mosi shale the same g function. Goeit/, et al. (it)7'l) found that 


= - H X 10 


ji)Ni;s,mi;lvii-i.i;: I'Ki-t.iMiNAHY moi)i;l srt)i)ii:s 


2i 


and hence 


Since they found 


b, = b^ Vl + (8 X iO- *p)*'. 


b. oc p- 


over the range p = 20 to p = 80, the corresponding p dependence for 
b^ should he correct eil hy the factor p“- ", or 


fc,ocp-«« 


so that 


a = 1.78, 
b = -0.01. 

As a check on this, we determined the power law depenilence of b^ on p 
directly and ohtuined values for a, A, and b of 1.75, 1.0 X 10', and +0.02 
respectively for the range p = 30 to 80. Coinhining our results with those 
of Goeilz, et ul. (197-1) we find that /j is given by 

/, = - eosh z/D + 15) (:W) S p S 80), 

where b has been assumed tx|Oal to zero, and, based upon one well- 
defined current dip, Du has been set erpial to 1 in units of Jovian radii, 
'fhe total F function representing the axi-syinmetrie portion of the field 
is then 


1.0 X 10^ 


(p» + z*)"* 


(log cosh z/D + 15). 


A plot of F is shown in Figure 2. Applying the same conditions as for 
the spherical model, the above model predicts tliat the last closed field 
line will cross the magnetic rxpiator at p^. IflO. Using 0=15 aiul.A = 
7.5 X 10 ', Goertz, ct al. (107-1) obtain p,, s 150. 

The Currents 


3'be current configuration in the magnetosphere can be obtained simply 
from Ampere’s law. Using the field expressions derived from the several 
fi functions one can obtain the configuration of the intense current 
sheet that exists at the magnetic erjuator as well as the vohiim- currents. 
The tf> components of the internal magnetospheric current svsleiu is found 
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I'lfMiif 2 . Mt'iiilmiial [il.iiic ir|iirM-iitaliiiii iil lilt' iiiii(;in'lii‘'|iln'riu lirlil its ilfvflii(ifil 
111 .1 t-vliiiiii i< III riMiiiliii,ili* rr|iirM'iil<itMiii. Tlir siilitl mill il.islii'tl iirriiws ri'|ni*M‘iil 
llir .ivri.i^i- lii lil iliirriiiiii iiiriiMirrtI iil si-vn.il |iiiliils .iliiiiH llii- iiiillHiiiiitl liiiji-t lllrll■^ 
III rioiiri I 10 mill II rrspi I livriy. Till- \li.iilrtl liiill wiillli |iiiilMiii iii-ar llii- iiiii)>iifti<' 
|'||II||I||| ii'|iirM'iils II |IIIII|IIII III llir 11111 . 111111.11 riiiii'iit slii'i'l riiiili^iitiiliiiii .ismiiiii-iI ill 
tins iiiiiilrl, anil a llisl iiiilri a|i|i|iixiiii.iliini In llii' iiia)>iirlii|iansi' iMiiiinl.iry i\ alxii 
Iiiiln .Ill'll. Ill y.iiiiix Iivi'i wliii li llir Inin liuii ix ri lialili- mi- niilii ali'il in llii- li'Xt. 

I hr I nivrs Ir.iv r ihr I' “ 2 S|ilirir ill rijiially X|iarril mi^nliir Inlrrviilx. 


Ill I'lirli rasi' In niiiMSI ill a slii'i'l ( iini'iil li'iiii plus a voliiiiic ciirn'iit 
Iriiii. w Ill'll’ llii' slii-i't li'iiii liil llir spliciical inuili'l is 

/\siiiW ^ 

' i lls «,i i iisflii 

.mil liii llii- 1'\ liiiili ii'al iiiiiilrl 




A 

Mil /MV'" 


si'fli- zH). 
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Allhiingli llii: voliinir torin.s iiri.* negligililr near the magiii'lic c(|iiatiir, 
they domiiiati' iwar tlir inagnclic polar axi.s. This is cli'aily an artifarl 
of rai-li iiituli'l wliii'ii will disappear when proper aeeoiiiil of ihe iiiagiii'lit- 
paiise cnrri'nts and orihe inner eiiloir radios of ihe eniri'iil shei'l or ilise 
are inehiih'il. 

Mallu'inalii.'ully (i'rniinating ihe field at the inagiu'lopaiiM' allows one 
li> .solve for (he inagnelopanse currents, and (he corresponiling honnilary 
iiclil ilirection can he ilelennined for the several niiiili ls anil coinpari'il 
with Ihe ilata. For example, a IniiL'Ilon, t{plpn), which can li'iininale thi' 
field arhilrarily ahriiptly is 

.,11 * + hmh r/(l - plp„) 

Uplpo) =* 2 * 

so (hat the terininateil fii'ld, // ', is given hy 

B' = tipMB. 

Ilerep,) is the railial ilislanee to (hi* magnetopause (as a liisl ap|>roxima- 
tion we assume the magni'topanse hmindary to he sphiTical) anil 1 / ri'- 
lules to (he thickiii'ss id' (he hoiindary. In principle, li slionld he ihi' total 
lielil. However, we still neglect the component of ihi' pi'ilnrliation 
lii'hl. 

'I'he ahove runction (ermiiiiites the preeeiling a/.imnthal iiirri'iits at 
p = po anil in adililion proviiles (he magnetopause cniTi'iits, i.e., for tlu' 
spherical model 

hi 

Ln ^Pii‘hPn 


and for thecylimlrical model (herep = Vx^ + y^, normali/.i'il) 



-il _ 
iPiAPn 




zli^-pU, 1 

Vp^ i xJ J' 


We liiid (hat li, for (he sphi'rical moilel is positive at all vahii's of 0, so 
that the corresponding imigni'lopanse current is cloi-kwise, as vii'weil 
from the magnetic puli', at all latitudes. Hence, just prior to the magneto- 
puiisu lionndury (he predicted field direction is sonthward, as is ohseiAcd 
hy hotli rioneers 10 and II (Smith, et al., 107-1; 107.’>). On the other hand, 
(he hrueketed term eoiitaiiied in the magnetopanse expression for (he 
cylindrical model heeomes negative at magnetic latitudes greater than 
ahont 20°, so (hat the ilirection of (he magnetopause enrn-nt How ri'- 
verses from a i-loi-kwise direction at lower latitudes to a eomiteri loek- 
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wise (lir<-(.-t ion at lii^licr laliliiilcs. Tlic conrspoiidiiig ilirld just insido 
llu; Ixuindaiy is predicted to point uinlhwiiiil at latitudes greater tliaii 
20° and sootliward at lower latitodi's. Snell a prediction ap|>ears to lie 
in disagreement witli tlie data, altliongli it is inteiesting to note tliat prior 
to tlie (iiitluaind inagnelopanse crossings liy I'ioneer 1 1 tliere were inter- 
vals approaeliing 8 hours during which the inagnetos|iheric Helds pointed 
noith ol radial liy roughly 20°. However, these are likely transient lea- 
tures related to a (possibly) tiorlhward component of the solar wind How 
vidocity. 

Discussion 

Sinci- the runetions jilotted in I'igures 1 and 2 wer«- th-rived from the 
rioneer 10 outbound data, they were devcIopcHj from data taken within 
■iliout 20" of the magnetic eipiator and over the radial range 20 = p = 80, 
and (pialitalively represent the inagnetospheric field conrignration in 
meridional planes that lie la-ar the dawn terminator. However, one 
notes that the models also ipialitatively fit the I’ioneer 1 1 data ipiite well, 
which e.\t<-nds the latitude range of tin' functions to peihaps ‘10° (the 
I’ioneer If inbound data is ipialitatively veiy similar to the I’ioneer 10 
outbound data) anil to .ibout 10° sunwaril of the dawn mei'itlian. As is 
evident from the figures, both models should be considered unreliable at 
latiliides gic.ilei than about ‘l.'>°. 

A basic dillcrcnce between till' two models is the fact that one is for a 
constant angular width current sheet (the spherical coordinate model, 
I'igme I ) while the other is lor a constant thiekness curient sheet (the 
cvlindiic.il coordinate model, l''ignre 2). I.ikely the actual ease lies .sonie- 
wlicic between these two current sheet configurations. Until permit the 
CIO I cut sheet to exist to the center of the planet, although it must b(‘ i-ut 
oil at Slime mliel radius p~ 2, since one would not expect the sheet to 
exist wilhiii the centi ilugal-graxitational balance distance of several radii. 

.'\uollier basic diflerence involves the dircetioii of flow of the magneto- 
p.mse ciiiients and the eorresponding direction of the magnelO|)aitse 
field. I'ledietions b.ised upon the spherical model are more consistent 
with the measmenienis. 

I here ,iie .ilso a number of factors regarding tlie constants derived for 
the models (hat should be mentioned. I'or example, in the case of the 
c\ lindiic.il ( ooidmale model, the value ol r in /, is deteriiiiiieil from 
/», x'cisns p at constant zID, but this nspiiri-s a knowledge of D. 'Hie 
ex ahi,ition ol tin- constants ii, ,A, and h depends critiealiv upon the aeeii- 
i,ite deieiniinalion ol the actual jiower law depeiidenee ol /) on p. llie 
I'oiienl sheet hall width is one ol the most miceitain patamelers, and its 
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dependenee nponp is piirtieiihirly difliciilt to determine ilireclly from the 
plots of li versus p. 'Hie cunstiint.s contained in the expri-ssion for J\ 
are self-consistent, anil the model appears to establish ihep imlependence 
of D. Similiir comments can be made regarding the i-onstant angular 
width model as well, tinforlimalely, a brief stuily of the variation of the 
widths of the lielil dips has not shed inncli light on this ta iieial [xiinl, ex- 
cept that the data lend to favor the consttml angular width model. 

In the study by Goerl/., et til. (1974), the ilelennination of n and D (their 
/>,)) from a |>lot of h, instead of causes the resulting function that is to 
represent the shape of the Held in meridional planes to be a mixture of 
f functions requiring dilferent g functions. On the other haiiil, in our 
delerminalion of these constants for the cylindrical model, we have used 
oidy f funelions that have the same g function, neeause the S|)iraling of 
the held was not excessive, the disagreement with the results of Goert/., 
el al. (1974) is not great, and a comparison of the plots of the held lines 
shows them to be quite similar. 

Any interpretation regarding the value of p,. (where Hi ~ 0) that is 
derived from the models should be viewed with caution since the model 
fils the ilala only onl to alMtuI 80 or 8.5 Hj, and hence these cutoff radii 
should be considered as possible artifacts of the models. An artifact of 
this kind is meaningless because the magnelopaiise currenls have been 
neglected in the derivation of the / functions. As noted earlier, it is 
tempting to assume that held lines leaving the planet at higher magnetic 
latitudes than those related to p,. are open field lines and that they merge 
with the inleqrhmelary field. But the data show that the field lines are 
.southward at the magnetopause, suggesting that they are closed by the 
magnetopause currents. In the sense of field lines and particle trapping, 
the.se lines clearly will not have trapped particles on them. The last 
closed field line which could contain trapped particles should be the one 
which crosses the eipiator just prior to the magnetopause boundary. 

'fhe particles in the intensi* equatorial current .sheet likely result from 
plasma How line to the combined action of a jovian "polar wind," mneh 
like that postulated for Karth (Banks and Holzer, IfXib), plus the strong 
centrifugal force caused by the large size and rapid rotation of the mag- 
netosphere. The balance of [rressiircs at the magnetopanse likely must 
include that exerted by a radial How of polar wind ions moving parallel 
to the e.s.sentially radial held lines in the magnetosphere. I’eihaps .such 
a plasma How also provides a significant stabili/ing inlhieuec for the 
large-scale magnetosphere conhguration reported here, since one would 
otherwise expect the .solar wind to blow the high latitude held lines 
back into the tail because of the relatively weak magnetic pressure 
exerted at the magnetopause (.Smith, et al., 1971). A study of the I’ioneer 
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11 oiitl>oiiiuI (lain will provide some impoilant infoniiation in tliis regard, 
ultlioii^li iiiiidi liiglier laliinde data arecleaiiy needed. 

Fiiillter stndi(!S of die iiiagnelo.splierc will likely reipiire (lie use of 
peiliirlialioii leeliniipies (Sleni, |!Xi7) in order to olilain more well- 
liehaved riinctions at high lalitnde.s and to allow for u non.sphirrieal 
inagiu'lopanse hoimdai'y. Cither sliidie.s heing condneted al the [iresent 
tiiiie will merge models developed for the range I ^ p S Cl with die 
magn(;los|)heiie models reported here. In this regard, magiurtospheric 
.stiidies estahlish reasonahle estimate.s of the magneto.spheriu eurrent 
systems and detailed attempts at merging the two piograms will estahlish, 
among other things, the inner enloif radius of the current sheet. 
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